Unravelling the glue and the competing order in superconducting cuprates 
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We present Raman scattering experiments in La2- x Sr x Cu04 single crystals at various doping 
levels x and compare the results with theoretical predictions obtained assuming an interaction 
mediated by spin and charge fluctuations. The light-scattering selection rules allow us to disentangle 
their respective contributions. We find that the glue spectral function is spin-dominated at low 
doping while the contribution of charge fluctuations becomes dominant around optimal doping. 
This indicates that the fluctuations of a nearly ordered state with coexisting spin and charge order 
support the superconducting pairing. 

PACS numbers: 74.72.-h, 78.30.-j, 74.20.Mn, 71.45.Lr 
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I. INTRODUCTION 



Several systems, ranging from heavy fermions to 
manganites, from high-temperature superconducting 
cuprates to ruthenates, display anomalies in the metallic 
phase, with the violation of the well established paradigm 
of Fermi-liquid theory. As a matter of fact all these sys- 
tems are characterized by strong electron-electron cor- 
relations and display similar phase diagrams, where an 
ordered phase is adjacent to the metallic state, possi- 
bly with the occurrence of superconductivity. This raises 
the question about the effective electron-electron inter- 
actions determining the metallic anomalies and possi- 
bly (high-temperature) superconductivity. In this re- 
spect cuprates are a paradigmatic example: Around op- 
timal doping (where the superconducting critical tem- 
perature T c is highest) these systems display anomalies 
in spectroscopic, transport, and thermodynamic proper- 
ties. These anomalies should stem from the same ef- 
fective interactions which also provide pairing. Under- 
standing their nature is the hotly debated "glue issue". 
For some people 1 - - — the glue consists of strong, essen- 
tially instantaneous, electron-electron interactions, aris- 
ing from the doped Mott insulator character of these 
systems. According to a second point of view, the low- 
energy effective interaction is due to a retarded bosonic- 
like glue 4 -^. Besides the standard phonons, spin waves& 
are the most quoted candidates, as reminiscent of the 
antiferromagnetism suppressed by doping. Quite natu- 
rally, the bosonic excitations might also be related to 
an elusive electronic order, which competes with super- 
conductivity and should occur in underdoped cuprates 
(i.e. at doping less than optimal) either as a long-range 
ordered phase or in the form of local/dynamical short- 
range fluctuations^—. The proximity to the correspond- 
ing "critical region" naturally brings along abundant 
critical fluctuations and leads to strongly momentum-, 
temperature- and doping-dependent bosonic effective in- 
teractions. Within this second point of view the iden- 
tification of the glue also sheds light on the underlying 



competing instability. 

Here we show that exploiting the specific properties of 
Raman spectra in the different symmetry channels, fluc- 
tuations with different characteristic wavevectors can be 
identified. This is a quite general result for effective in- 
teractions that are strongly peaked at a finite wavevector 
(see Appendix for technical details) and it might render 
Raman spectroscopy a particularly powerful tool to de- 
tect and study elusive orders occurring at finite wavevec- 
tors. More specifically, here we start from a phenomenol- 
ogy for spin and charge fluctuations, predict the Raman 
spectra in the two Bi g and £?2 ff symmetry channels and 
compare our theoretical results with experimental results 
on La2- x Sr x Cu04 (LSCO) samples at various doping lev- 
els x and temperatures T. In this way we disentangle the 
contribution of spin and charge fluctuations, which can 
be identified as the relevant bosonic modes in these sys- 
tems. This identification, not only clarifies the glue issue, 
but also strongly indicates a nearly spin/charge-ordered 
state as the competing phase in underdoped cuprates, 
which evolves from a dominating spin order at low dop- 
ing, near the Mott insulator, to a charge-dominated order 
on the overdoped side. 



II. EXPERIMENTAL DETAILS 

The samples were grown using the traveling solvent 
floating zone (TSFZ) method. The characterization is 
presented in Table fl] 



The lowest doping was x = 0.15 close to the T c maxi- 
mum where Landau-Fermi liquid theory should still be 
applicable. The overdoped sample with x — 0.26 had no 
indication of superconductivity above 2K. All samples 
were post-annealed in pure oxygen at 1 bar to improve 
the crystal quality. The spectra were taken on polished 
surfaces. In the case of Lai^Sro^CuCU the results were 
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TABLE I: List of samples. Samples labeled with a have been 
prepared by M. Lambacher and A. Erb (WMI Garching}^, b 
by Seiki Komiya and Yoichi Ando (CRIEPI, Tokyo and Osaka 
University) and c by N. Kikugawa and T. Fujita (Hiroshima 
and Tokyo). The transition temperatures were measured ei- 
ther resistively or via magnetometry or via the non-linear ac 
response. 



sample 
Lai.85Sro.i5Cu04 
Lai.8 3 Sro.i7CuC>4 
Lai.8oSr .2oCu04 
Lai.7 5 Sr .25CuO4 
Lai.74Sro.26CuC>4 



doping x T c (K) AT C (K) comment 



0.15 
0.17 
0.20 
0.25 
0.26 



38 
39 
24 
12 




O2 annealed a 
O2 annealed b 
O2 annealed a 
O2 annealed a 
O2 annealed c 



compared to those from a cleaved surface and were found 
to be identical to within the experimental error—. 

The Raman experiments were performed with stan- 
dard equipment. For excitation an Ar ion laser operated 
at 458 nm was used. The angle of incidence of the incom- 
ing photons was 66°. The polarization state outside was 
prepared in a way that photons inside the sample had 
the desired state. Scattered light of a selected polariza- 
tion was collected along the surface normal and focused 
on the entrance slit of the spectrometer. The energy- 
selected photons were detected with a liquid nitrogen 
cooled CCD detector. All spectra are corrected for the 
sensitivity of the complete setup. We generally show the 
imaginary part of response functions, x" at pure sym- 
metries fj, = B\g, i?2 9 which are related to the measured 
cross section as 



(TOjj, 



n( W ,T)} X "(q = 0,w,T) (1) 



with cr M the photon cross section in symmetry fi, Q the 
solid angle accepted by the collection optics uj — uii — 
w s the energy transferred to the system, A a constant 
absorbing all factors to convert into a cross section 
and n(uj, T) = [cxp(w/T) — 1] _1 the Bose-Einstein factor. 
To simplify the notation we drop hereafter the index /z, 
since the formal expressions in the following are the same 
for all symmetry channels. 



III. RAMAN SPECTRA 

Electronic Raman scattering is a bulk (nearly surface- 
insensitive) probe and it measures a response function 
x"(w,T) analogous to that of the optical conductivity 15 . 
However, while the latter yields an average over the entire 
Brillouin zone, light scattering projects different parts of 
the Brillouin zone for different polarizations of the incom- 
ing and outgoing photons^, thereby introducing specific 
form factors in x"(u,T). We consider Raman spectra in 
the B\g and Big channels obtained from the systematic 
analysis of several LSCO samples at different doping lev- 
els x ranging from the nearly optimal up to the strongly 



overdoped region (see Fig. [T]). Generally we observe that 
the overall shape of the spectra in B\ g and Bi g channels 
is quite different up to 3000-4000 cm -1 , while at higher 
frequency the two spectra have similar shape and inten- 
sity. At frequencies up to 2000 cm" 1 , the Bi g spectra can 
be schematized by a hump appearing as a shoulder of a 
slightly larger hump peaked at to ~ 4000 cm -1 . On the 
other hand, Big spectra display a plateau up to about 
1000 cm -1 . At higher frequencies the absorption rises 
and a hump nearly twice as large, similar to the one found 
in Big, appears at about 4000 cm -1 . The temperature 
dependence of the spectra is altogether weak (see Fig. 
[2^), but it may even be opposite in the two channels in an 
intermediate frequency range. In particular the tempera- 
ture dependence of the Bi g spectra corresponds to those 
of the optical conductivity for x > 0.05. In Bi g symme- 
try at low-intermediate dopings there is a strong extra 
contribution from fluctuations at low frequenc y 43 ! 44 ] 49 . 
Therefore we confine the analysis to x > 0.15. 

A first remark is now in order: Although the Raman 
form factors select different regions of the Brillouin zone, 
and therefore explore different electronic structures, it 
is rather unlikely that the changes in doping and tem- 
perature in the two channels [see the low-energy part of 
Figs. 2(a) and the overall shapes of the responses in Figs. 
l(b-e)] can simply be attributed to diverse evolutions of 
the underlying fermionic band structures. Also a scatter- 
ing mechanisms with zero characteristic momentum^, 
would act in the same way in the two channels hardly ex- 
plaining their different behaviours. Guided by a wealth 
of experimental evidences for a prominent tendency of 
LSCO cuprates to form charge-ordered states (stripes or 
checkerboard) , we show below that different shapes (and 
different temperature and doping evolutions) of the spec- 
tra in the two channels naturally arise from quasiparti- 
cles (QPs) coupled to distinct charge and spin fluctua- 
tions, each one having its own dynamics and character- 
istic wavevector. 

IV. THE "GLUE" COLLECTIVE MODES 



Here we approach the identification of the glue extract- 
ing its properties from the Raman data. We start from 
a form of damped glue boson collective modes (CMs), 
which is customary for spin fluctuations^^ 



#A(q,w) 



1 



m\ + ^A(q- qAj 



(2) 



where A = c, s refers to charge or spin CMs substantially 
peaked at characteristic wavevectors q s ps (71", tt) and 
q c fa (±7r/2, 0), (0, ±7r/2) respectively. Here the mass 
m\ is the minimum energy required to excite the CM 
and v\ is a fermion scale setting the CM momentum dis- 
persion. This dispersion is limited by an energy cutoff 
A^. The dimensionless quantities m\/ K\ are the inverse 
square correlation lengths (in units of the lattice spac- 
ing), which measure the typical size of ordered domains. 
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FIG. 1: (a) Experimental Raman response ^u(w,T) of LSCO at x — 0.20 on a broad frequency range up to 8000 cm 1 at 
intermediate T = 104 K in the B\ g channel (green circles) and in the Bi g channel (red triangles), (b-e) Raman response 
X^(w,T) of LSCO at x = 0.15,0.17,0.20,0.25, and 0.26 on a reduced frequency range up to 4000 cm" 1 . Data in the B lg 
channel are given by green circles, while the spectra in the B2 g symmetry are given by red triangles. The various spectra are 
taken at similar temperatures around T 100K. The results of the fits obtained within the "nearly-critical" charge/spin theory 
(see text) are given by the black solid lines. 




FIG. 2: (a) Raman response T) of LSCO at x = 0.20 in the Bi g (circles) and B^ g (triangles) at three different temperatures 

above T c and the related fits (solid curves) obtained within the nearly critical scheme described in the text: T — 56 K (blue 
circles and curve), T = 104 K (green circles and curve), T = 201 K (red circles and curve), (b) Temperature dependence of 
the inverse square correlation length £,7 2 of the charge modes at different dopings (x = 0.15 blue, x = 0.17 red, x — 20 green, 
x = 0.25 black), (c) Same as in (b) for the inverse square correlation length £J 2 of the spin collective modes, (d) Temperature 
dependence of the low-frequency typical energy scale uii = V mQ. of the charge and spin collective modes. 
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FIG. 3: Hot spot structure for both charge (with q c ) and spin 
(with q s ) scattering. In panels (a) and (b) the sign structures 
of the form factors for the Bi g and B2 g channels, 7i s (k) = 
cos(fc^) — cos(ky) and 72 9 (k) = sin^) sin(k y ) (we use a unit 
lattice spacing), respectively, are reported. 



The iuj term establishes the low-energy diffusive charac- 
ter of these fluctuations due to decay into particle-hole 
pairs, whereas above the scale set by £l\ the CM has a 
more propagating character. 

The CMs mediate strongly momentum-dependent 
scattering between the QPs. At low energy this scat- 
tering mechanism is more effective for QPs on the Fermi 
surface, which are connected by the characteristic CM 
momenta (the so-called "hot spots"). For the typical q s 
and q c in the cuprates the hot spots for spin and charge 
CMs occur in the same regions of the Fermi surface near 
the (±7t, 0) and (0, ±7r) points of the Brillouin zone (see 
Fig. El). 

The Raman response function can be expressed in 
terms of a memory function as 



M{lo) 



(3) 



Here xo 1S the purely real low-frequency Raman response 
function in the absence of any scattering process. As in 
standard approache d 18 ' 19 , our perturbative memory func- 
tion calculation is limited to processes which involve the 
exchange of one CM (A = c or s). The specific form of 
the fluctuations ([2]) allows one to identify and analytically 
calculate the dominant contribution to M = M c +M s (for 
details see the Appendix) yielding for the imaginary part 



ImM x (uj) 



1 



dz [a 2 F x (z)] 



2ui coth ( — — ) 
\2TJ 



. z + uj\ , / z — u) 

[z + d) coth + (z — uj) coth 



2T 



2T 



4) 



which has the meaning of a frequency-dependent inverse 
scattering time. The crucial physical ingredient in Eq. 
(|3]) is the spectral distribution a 2 F appearing, e.g., in 
the textbook Eliashberg theory of superconductivity (the 
so-called "glue function" ) , which in our case reads 



a 2 F x {uj) = g\ 



arctan 



arctan 



m\fl\ — oj 2 



(5) 



The dimensionless parameter g\, which also absorbs v\, 
describes the coupling between CMs and QPs. This spec- 
trum is roughly contained in the range between u>\ ~ 
VOrn and ui2 ~ When fl is large, a 2 F(uj) has 

the broadly peaked shape of a glue function due to diffu- 
sive CMs. For small f2, i.e., for rather propagating CMs, 
a 2 F(oj) has a "flatfish" box-like shape. This latter form 
of the scattering spectrum has been phenomenologically 
assumed since the early days of high-T c superconductiv- 
ity by the marginal-Fermi-liquid theor y 20 ' 21 . We find it 
non trivial that the marginal Fermi liquid phenomenol- 
ogy in momentum-integrated quantities (tunnelling spec- 
tra, optics, Raman, etc.) is obtained here from a micro- 
scopic strongly momentum dependent scattering mecha- 
nism. When adapted to the calculation of optical con- 
ductivity, our formalism contains as limiting cases the 
marginal-Fermi liquid and the diffusive CM considered 
in Ref. and reproduces their results. 

Our first general result is that the different experimen- 
tal spectra in the B\ g and Bi g channels are obtained and 
correspond to a more diffusive and a more propagating 
CM respectively (an extended presentation of this generic 
result is given in the Appendix) . 

V. ANALYSIS OF THE COLLECTIVE MODE 
SPECTRA 

Despite the fact that both charge and spin CMs cou- 
ple to the QPs in the same hot-spot regions of the Bril- 
louin zone, symmetry arguments imply that their dif- 
ferent characteristic wavevectors gives rise to important 
cancellations in the memory function depending on the 
channel form factors. The cancellation occurs when the 
relevant scattering wavevectors connect regions where 
the form factors have the same sign (see Eq. (IA1|) in 
the Appendix). Conversely the largest contribution is 
obtained when the critical wavevectors connect regions 
with opposite sign of the form factors. Inspection of Fig. 
[3] shows that the spin CM dominates the B\ g spectra, 
while the charge CM dominates the spectra in the B^ g 
channel. These "selection rules" are effective at low fre- 
quency and gradually fail upon increasing the energy of 
the Raman scattering. This failure is effectively taken 
into account by a gradual frequency-dependent switch- 
ing on (from about 1000 cm -1 ) of the coupling between 
the QPs and the "forbidden" CM in each channel (for 
details see Appendix ). We do not attempt to fit the 
data above 4000 cm -1 because multiple scattering pro- 
cesses become important making the distinction between 
M c and M s meaningless. At this point the spectra even- 
tually become similar in shape [see Fig. HJa)]. 

Adjusting the CM parameters and their coupling to 
the QPs we fit the experimental spectra (see Fig. [TJ). 
Remarkably, the different shapes and doping evolutions 
can be tracked. For instance, the B\ g spectrum evolves 
from the rounded shape (characteristic of a strong scat- 
tering already at low energies) at x = 0.15 to the nar- 
row peak typical of a weakly scattered "Drude-like" be- 
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haviour at x = 0.26. From these systematic fits we ex- 
tract the corresponding evolution of the glue functions 
(see Fig. EJ and of the related CMs [see Fig. H>d)]. 
Quite naturally in the substantially doped regime we are 
considering, we find that the spin spectra (dashed lines) 
have a marked diffusive character (typical values of SI are 
in the range of fOOO — 3000 cm" 1 ). The overall weight 
W s of the glue functions due to spin exchange [see the 
dashed lines in Fig. 4(f)] is still rather large at x = 0.17, 
but then decreases and becomes very small in the most 
overdoped sample (x = 0.26) in agreement with inelastic 
neutron scattering data2&2i. Conversely, the shape of the 
charge-mediated glue is more propagating-like (roughly 
resembling the shape proper to a marginal-Fermi-liquid 
glue) because the typical SI ~ 200 cm -1 is much lower 
than in the spin case. Thus, the charge-mediated glue 
function is centered at typical phonon frequencies of the 
cupratcs. This agrees well with the idea that charge- 
ordering in cuprates is phonon-driver* 7 -, with a strong 
mixing between the electronic and the lattice degrees of 
freedom 25 . Charge CMs have a marked phononic char- 
acter at generic q and acquire a strong electronic compo- 
nent at q ~ q c , where this composite excitation has an 
anomalous softening with weight displacement down to 
low energies of order m c . This behaviour is reflected in 
the temperature dependence of the charge glue functions, 
which, while having their centre at about 500 cm , also 
have a substantial spectral density at low energy when 
temperature is low. 

The overall weight of the charge glue W c evolves in 
the opposite way with respect to the spin weight. It is 
relatively small at x — 0.15, increases reaching a maxi- 
mum at x = 0.20 — 0.25 and then slightly decreases at 
x = 0.26 [see Fig. 4(f)]. Of course this doping evolution 
should be reflected in the doping evolution of the super- 
conducting critical temperature. Owing to the strongly 
retarded character of the interaction, a full Eliashbcrg 
analysis of superconducting properties is required. 

Both charge and spin CMs have strongly temperature 
dependent inverse square correlation lenghts £~ 2 = m/A 
(see Figs. 2(b), (c)) which only affect the low-frequency 
part of the glue functions, but have a marked influ- 
ence on the Raman spectra over the whole frequency 
range. Without this temperature dependence (like, e.g., 
for standard phonon scattering), the spectra would be 
much more temperature dependent, due to the bosonic 
occupation factors in Eq. 

In general, we find that the observed temperature de- 
pendence of the charge inverse square correlation length 
£~ 2 is in qualitative agreement with an underlying criti- 
cal behaviour related to a competing charge-ordering in- 
stability. In particular, we notice that the linear tem- 
perature dependence at x = 0.20 extrapolates to zero at 
T — > indicating that this doping is quite close to the 
value of the quantum critical point usually reported to 
be at x = 0.191&2&). 

We also find that the spin correlation length is of the 
same order as the correlation of the charge ordering at 



low T, consistent with the idea that spin fluctuations can 
survive at high doping, far from the antiferromagnetic 
region, if sustained and "enslaved" by charge-depleted 
fluctuating regions. This connection is also apparent in 
Fig. [2jd), where a similar energy scale uj\ for charge 
and spin is reported. This shows that the "bulk" of the 
spectrum for the two collective degrees of freedom has 
essentially the same energy threshold. 

VI. OTHER EXPERIMENTS 

Coupling functions derived recently from the Raman 
spectra in one Bi22f 2 sample yield differences in the two 
symmetries, with a prominent spin scattering in the B\ g 
channel 27 . However, the identification of spin and charge 
modes and the assessment of their contribution to the 
two channels at low frequency has been made possible 
by the procedure described here. The present work may 
solve the long-standing puzzle in the cuprates regarding 
the relevance of the spin vs. phononic scattering mech- 
anism affecting, e.g., the electron dispersion and giving 
rise to the kinks in ARPES spectr a 6 ' 28 ! 29 : both spin and 
mixed charge-phonon modes contribute to the anomalies 
of the electronic dispersion because both modes share the 
same hot-spot region and have a substantial weight. It 
is also encouraging that recent ARPES 30 and STM21"^ 
experiments support the role of charge-order scattering 
in electronic spectra. Moreover, the relevance of charge- 
order excitations with a marked phononic character, can 
also provide a rationale for various isotopic effects de- 
tected along the years in cuprates^— 

The simultaneous presence of two CMs dressing the 
QPs, was also recently found in optical spectroscopy on 
Bi 2201 and Bi 2212 samples 3 ^. One could also sim- 
ilarly interpret the a 2 F recently obtained from ARPES 
spectr a i . These findings are in agreement with our 
analysis, because in optical conductivity and in ARPES 
the channel-dependent form factors are absent and both 
CMs contribute. As a result we simply need to add the 
spin and charge glue-functions to reproduce the exper- 
iments. Although we did not attempt a precise fitting 
and comparison with the experimental data, we explain 
why the optical scattering time (Jm M opt (w)) above 1000 
cm displays a very weak (if not absent) temperature 
dependence in underdoped samples, which gradually in- 
creases upon doping (see Fig. [SJa)-(c))^ 5 - _ — . This be- 
haviour simply arises from the temperature dependence 
of the CM masses (both charge and spin): upon increas- 
ing T, the masses grow producing a less marked increase 
of the scattering. This compensation is less effective upon 
increasing doping. 

While the above spectroscopies indirectly probe the 
CM through their effects of the QPs, other spectroscopies 
directly probing the charge or spin CMs provide pictures 
consistent with our findings. Concerning the spin modes, 
there is a correspondence between our spin glue functions 
and the spin spectral functions obtained from inelastic 
neutron scattering. Specifically we found that above 
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FIG. 4: (a)-(e) Glue functions due to the "allowed" modes in the fitted spectra of Fig. 1. The solid lines are the glue functions 
due to charge CM, while the dashed lines are the spin CM glue functions. Blue, green and red curves refer to low (T ~ 50 K), 
intermediate (T ~ 100 K), and high (T ~ 150 - 200 K) respectively, (a) x = 0.15, T =42 K ,95 K, and 195 K; (b) x = 0.17, 
T =54 K, 104 K, 254 K; (c) x = 0.20, T =56 K, 104 K, 201 K; (d) x = 0.25, T =53 K, 104 K, 201 K; (e) x = 0.26 T =46 
K, 102 K, 149 K; (f) Total weight of the glue functions mediated by charge W c (solid lines) and spin W s (dashed lines) as a 
function of doping x 



(b) x = 0.2C l^{c) x = 0.25 




2000 2000 2000 
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FIG. 5: Imaginary (a-c) and real (d-f) parts of the optical 
memory functions obtained summing the charge and the spin 
glue functions reported in Fig. [4] The correspondence be- 
tween colors and temperatures is the same as in Fig. U For 
details of the calculations see the Appendix 



x = 0.20 the spin excitation intensity rapidly decreases, 
by about a factor one half between x = 0.25 and x = 0.26. 
This corresponds well to the spin suppression observed in 
Refs. 23.24. Dynamical charge fluctuations are less ac- 
cessible. However, we interpreted previous results of Ra- 
man absorption in underdopcd LSCO samples at low fre- 
quencies (up to a few hundreds of wavelengths)^ in terms 
of direct excitations of two charge CMs^. Remarkably 
the values and the temperature dependencies of m c (T) 



obtained at such low frequencies and in much less doped 
samples, and from different excitation processes, are in 
close agreement with the behaviour of m c (T) found in the 
presently analyzed optimally and overdoped samples. 



VII. SUMMARY AND PERSPECTIVES 

In this paper we demonstrate that a retarded bosonic 
glue accounts well for the rich variety of shapes, and of 
doping and temperature dependencies of Raman spec- 
tra in LSCO. We show that, fully exploiting the inter- 
play between the sign of the Raman form factors and 
the peaked momentum structure of the fluctuations, the 
dominant absorption in the B± g and i?2g channels at low 
energy mostly arises from QPs scattering either with spin 
or charge CMs respectively. The different shapes of the 
Raman spectra in B\ g and Bi g symmetry is the result of 
the more diffusive or propagating character of the corre- 
sponding dominant CMs. The more propagating charge 
CM corresponds to a rather flatfish glue function thereby 
reproducing the phenomenology of marginal- Fermi-liquid 
theory in momentum integrated responses. 

Our main result is a complete description of doping 
and temperature evolution of the glue functions due to 
charge and spin CMs. The relative importance of the 
two scattering mechanisms switches from spin to charge 
by increasing doping. While the relevance of dynamical 
spin scattering was already inferred both from theoret- 
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ical approaches^ and phenomenological analyses^, the 
simultaneous presence of mixed phonon-charge-density 
fluctuations and their increasing relevance with doping 
is assessed here for the first time. Moreover, the rele- 
vance of charge fluctuations also indicates that a nearly 
charge-ordered phase competes with superconductivity 
in underdoped LSCO. Our analysis opens the way to the 
wide field of investigations based on a specific form of a 
retarded bosonic glue. 
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Appendix A: Memory function 

We outline here the main steps of the calculation of the 
memory function that appears in the Raman response 
equation (2). To simplify the notation we drop the index 
/x, since the formal expressions are the same for all sym- 
metry channels. Within a perturbative approach limited 
to the exchange of a single CM, the memory function 
reads 



M 



Xsv 
> . 

Xo 



where xo is the purely real low-frequency Raman re- 
sponse function in the absence of any scattering process, 
and xsv is the contribution due to CM selfenergy (S) and 
vertex (V) corrections (see Fig. [6]). Here we are neglect- 
ing impurity scattering as well as sources of scattering 
other than CMs. The fermion loops entering the expres- 
sion of XsvI^m), where Q,m is the external Matsubara 
frequency, are 

/j2i_ 
^ ry £0(k,e/-IW) 

x g(k,e e )g(k + q,e e + u n )g(k,e e ), 
S 2 (q,0J n ,n M ) = -2T~g 2 Y, J 7^7k6(k,e<+%) 

V(q,u> n ,n M ) = -2T~g 2 Y, J 7^j2 7k7k+ q 0(k, e t ) 

x g(k,e e + fl M )G(k + q 1 e e +uj n + ri M ) 
x g(k + q,e e + w„), 

where g(k,eg) = (ieg — Ck) _1 is the Matsubara Green's 
function of free fermions with band dispersion £k- To 
calculate xsv, we make now two main approximations. 
First, we notice that the CM propagator is peaked at 
a wavevector q = qA, whereas the fermion loops are 



smooth functions of the integrated wavevector q. Thus, 
in the integral over q that appears in the expression for 
Xsv, the fermion loops can be factorized and evaluated 
at q = qA • Once this factorization is introduced, the inte- 
gral over k in the fermion loops is dominated by QPs near 
the hot spots (HS)at the Fermi surface associated with 
q> • The two-dimensional integral over k is then suitably 
transformed into an integral over the QP dispersions at 
the two hot spots connected by q\, £hs and £' HS (this 
approximation corresponds to a linearization of the band 
dispersion around the hot spots, see Rcf. 6 ). This change 
of variables involves the Jacobian of the transformation J 
evaluated at the hot spots, giving rise to a dimensionless 
prefactor J = Jg 2 . The subsequent sum over the fermion 
Matsubara frequency ee may be carried out explicitly. 

In the forthcoming discussion, it is crucial to ob- 
serve that the selfenergy contribution contains two Ra- 
man form factors evaluated at the same hot spot, "f HS ', 
whereas one of the two form factors appearing in the 
vertex contribution, "f HS , is evaluated at a different hot 
spot, connected to the former by q\. 

After evaluating the fermion loops, within the approx- 
imations outlined above, we obtain the contribution to 
the Raman response 



Xsv(M 



M. 



= Ty^T sv {Lo n ,£l M ) 



(2ny 



:V(q, w n ), 



where 2?(q, uj n 



-\m\ + v\(q - q\) 2 



LJ n /Cl\]~ is the Matsubara CM Green's function 
and the leading contribution of the fermion loops 
is fully included in the function Tsvi^ni &m) = 
[Si{q,w n ,£l M ) + S 2 {q, u n ,Q, M ) + V(q, ui n , ^«)] q=qx - 

By explicit evaluation, following the procedure out- 
lined above, we find 



r sv (w n ,n M ) = 47rj 7 ffS ( 7 
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.MS 



in 



M 



fi 2 



for — ]Oat! < w « < I^a/I, and r S v(w„,n M ) = else- 
where. This expression enforces the selection rule quoted 
in the text: the main contribution to the fermion loops 
is non vanishing only if the Raman form factor 7k has a 
different sign on the two hot spots connected by qA- For 
the characteristic wavevectors of spin and charge CMs, 
this occurs in the B\ g and B 2g symmetry channel, re- 
spectively. 

After a shift of the origin of momentum space to qA, 
the integral over q becomes straightforward. The sum 
over the Matsubara frequencies u„ is then transformed 
into an integral over the real axis by standard techniques, 
exploiting the fact that the boson frequencies ui n are the 
simple poles of coth(z/2T). After causal continuation 
of the external Matsubara frequency £Im to the real fre- 
quency oj, we finally obtain the contribution to the Ra- 
man response function due to the QP scattering with 
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FIG. 6: Self-energy (S) and vertex (V) diagrams corrections 
to the Raman response function with the exchange of a single 
collective mode. The black dots represent the 7k form fac- 
tor, while the shaded dot represents 7k+ q - The dotted lines 
represent collective mode Green's functions 
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where the symbol P indicates that the principal part of 
the integral must be considered whenever a simple pole 
of the integrand is met on the real axis. Although this 
is not immediately evident, by changing the integration 
variable to — z and summing one half of each of the two 
expressions it may be easily shown that the imaginary 
part of the resulting memory function is an odd function 
of the external frequency w, as required by causality. This 
explicit antisimmctrization with respect to uj yields the 
results for ImM and a 2 F quoted in the text, with the 
dimensionless effective coupling 



9\ = 



JjHS ^HS _ 7 HS'"j 



In the presence of multiple hot spots and/or multiple val- 
ues of the characteristic wavevectors qA (as implied, e.g., 
by the symmetry of the lattice), the overall dimensionless 
coupling is intended to be a sum over all the hot spots 
and qA's. 



Appendix B: Shape of the Raman response 

The connection between the overall shape of the Ra- 
man response and the properties of the CMs is fully en- 
tailed in the form of the CM glue function a 2 F(uj). In the 



FIG. 7: (a) Glue functions for two CMs having the same 
dynamical range between ui\ = 330 cm -1 and L02 = 1500 
cm -1 , but different diffusive-propagating crossover scale: ft — 
200 cm -1 for the more propagating mode (solid red curve), 
fl — 1000 cm -1 for the more diffusive mode (dashed black 
curve), (b) Imaginary part of the memory function (scatter- 
ing time) corresponding to the two glue functions in (a), (c) 
Raman response obtained for the two glue functions in (a) 
with g 2 = 0.5. 



forthcoming discussion it is important to observe that the 
overall behavior of the Raman response is qualitatively 
understood by inspecting the behavior of the frequency 
dependent inverse scattering time 1/t(cj) = ImM(uj), 
since the mass corrections associated to ReM{u) do not 
introduce significant features. 

A nearly constant 1/r produces a Drude-like peak 
at frequency u) ~ 1/t. A nearly linear 1/r produces 
instead a flatfish (i.e., marginal- Fermi-liquid-like) re- 
sponse. Now, we see that a diffusive CM (large O) is char- 
acterized by a broadly peaked glue function [see black 
dashed line in Fig. 0a)] . This, in turn, yields a 1/r which 
smoothly interpolates between the low-frequency value 
(l/r ) and the high-frequency value (1/too > l/r ) [see 
black dashed line in Fig. EJb)]. The Raman spectrum is 
correspondingly characterized by a hump+hump overall 
shape which results from the overlap of two Drude peaks 
roughly centered at l/r and 1/tqo [see black dashed line 
in Fig. 0c)]. This corresponds to the B\ g spectra ex- 
perimentally observed in LSCO. If the CM has a marked 
propagating character (small the corresponding glue 
function has instead the shape of a rounded box function 
[see red solid line in Fig. 0a)]. This, in turn, implies 
a broad linear regime in 1 /r and a wider separation be- 
tween 1/ro and 1/roo. The corresponding Raman spec- 
tra has the step+hump form [see red solid line in Fig. 
0c)], similar to the Pig spectra experimentally observed 
in LSCO. 
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Appendix C: Relaxing the selection rules 

The channel-dependent selection rule for the Raman 
spectra was obtained under two main assumptions: (1) 
the CMs are peaked at finite wavevectors and (2) the 
factorized fermion loops are dominated by QPs at the 
hot spots. The second assumption is asymptotically valid 
at low frequency, but is gradually violated at higher fre- 
quency. This means that at sufficiently large u> both CMs 
contribute to Raman spectra, regardless of the symmetry. 
This is expected on theoretical grounds and indeed ob- 
served experimentally: the B\ g and Bi g spectra of LSCO 
have similar shapes at high frequency, indicating that 
channel-dependent selection rules are not at work there. 

However, the analytical evaluation of the fermion loops 
beyond the hot-spot contribution is not viable and the 
numerical evaluation is both demanding and not trans- 
parent. Thus, in our calculation we adopted an alterna- 
tive phenomenological description, writing the memory 
function as M(uj) = M a (uj) + s(w)M/(w), where M a (w) 
is the memory function associated with the CM allowed 



by the selection rule and Mf(u>) is the memory function 
associated with the CM forbidden by the selection rule. 
The switch function 



«(w) = 



1 



1 + tanh 



gradually includes the contribution of the forbidden CM 
at a frequency cu s=s u>f. The rapidity of the inclusion 
is controlled by the parameter fij. Typical values in 
our fits are ujj <~ 10 3 cm -1 and /?/ <~ 5 x 10 2 cm -1 , 
roughly common to both charge and spin modes. We 
point out that our phenomenological procedure is not 
fully causal (i.e., ReM and ImM are no longer related 
by Kramers-Kronig transformations). However, since a 
constant s(u) would restore causality, a sufficiently grad- 
ual inclusion of the forbidden CM does not introduce 
significant drawbacks. Below the frequency w/, the spec- 
tra in the B\ g and Bi g symmetries are markedly differ- 
ent, whereas above w/ they gradually become similar in 
shape. 
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